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ABSTRACT 

Cold optically thick accretion disks with hot coronae and radial advection have 
been investigated. Within the framework of a-viscosity models, we assume that all the 
mass accretion and angular momentum transport take place in the cold disk, but that 
a fraction of the gravitational energy released is dissipated in the corona. Both the 
coronal energy dissipation and the advection heat transport have a stabilization effect 
on the thermal and viscous instabilities of the disk. If that more than ~ 95 percent 
of the total power is dissipated in the corona, then the locally unstable behavior of 
the disk is restricted to a relatively narrow spatial region and is found to lie in a 
small range of mass accretion rates. The global temporal variability of the disk can 
be very mild or may disappear, and which may be applicable to the low-frequency 
(~ 0.04 Hz) quasi-periodic oscillations observed in black hole candidates Cyg X-l and 
GRO J0422+32. 

Subject headings: accretion, accretion disks — black hole physics — instabilities 

1. INTRODUCTION 

The spectral and temporal behavior of galactic black hole candidates (BHCs) and active 
galactic nuclei (AGNs) contains valuable information about the underlying physics of the accretion 
process. BHCs such as Cyg X-l have very distinct characteristics. Cyg X-l exhibits low and high 
states as defined by the X-ray flux in 1-10 keV band. The luminosity difference between these 
states is approximately an order of magnitude. In the high state, the spectrum consists of two 
components, namely, a relatively stable soft blackbody component and a highly variable hard 
power-law component. In the low state, the intensity shows rapid chaotic temporal fluctuations 
and the spectrum is characterized by a power-law component over an energy range from 1 keV 
to more than 100 keV. For a detailed review see Liang & Nolan (1984) and Tanaka & Lewin 
(1993). Characteristic spectral properties of AGNs include the so-called UV bump blackbody 
radiation, the power-law spectra of X-rays and the presence of Fe Ka lines (see, e.g., Pounds et 
al 1990; Matsuoka et al 1990). It is now well accepted that the soft X-ray black body spectral 
component from BHCs and the UV black body radiation from AGNs are due to a cold optically 
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thick accretion disk. On the other hand, the power-law spectra from these sources reflect the 
presence of optically thin hot matter. 

There are two general scenarios proposed for the optically thin hot matter, namely, an 
entirely optically thin hot accretion disk or a hot corona plus an optically thick cold disk. Early 
work was published in mid-1970s, such as Shapiro, Lightman, & Eardley (1976) for disk models, 
and Liang & Price (1977) for corona models. The structure of coronae has been studied with 
different heating mechanisms for the coronal gas. For example, Compton scattering of X-rays 
from the boundary layer or the innermost regions of an accretion disk may lead to the formation 
of a hot corona above the disk and a wind ( Begelman, McKee, & Shield 1983). Recently, Haardt 
& Maraschi (1991) proposed a two-phase thermal disk-corona model for Seyfert galaxies. In this 
model, most of the thermal soft photons emitted from the optically thick cold disk escape from the 
system by passing through the optically thin hot corona. A small fraction of them are Compton 
upscattered in the corona by the high-temperature thermal electrons, and which produces the 
hard power-law spectral radiation. About half of this radiation is directed toward the disk (i.e., 
irradiation) and is reprocessed to emerge as blackbody radiation; while the remaining half escapes 
from the system and is observed. The corona is characterized by its Thomoson scattering optical 
depth, r es , and the fraction of gravitational energy dissipated there, tj. This energy is balanced by 
the cooling of inverse Compton scattering. In order to obtain a power-law spectrum with an index 
of ~ 1 observed from Seyfert galaxies, it is required that tj ~ 1. Along this line, the structure of 
disk-corona system has been studied, for example, by Nakamura & Osaki (1993), Kusunose & 
Mineshige (1994) and Svensson & Zdziarski (1994) with different assumptions regarding such as 
the mass accretion and angular momentum transport taking place within the cold disk and corona, 
and the irradiation heating on the disk. For example, Svensson & Zdziarski (1994) assumed that 
mass accretion and angular momentum transport take place only in the disk, and no irradiation 
heating on the disk is included. 

Optically thick accretion disks have been widely studied. It is well known that accretion 
disks may suffer thermal- viscous instabilities when radiation pressure is important. In a standard 
a-model (Shakura & Sunyaev 1973), the viscous stress, r, is proportional to the total pressure, p, 
i.e., t = —ap, where a is a constant. Disks are thermally and viscously unstable if the pressure is 
radiation dominated (Lightman & Eardley 1974; Shakura & Sunyaev 1976). The thermal- viscous 
instabilities may be stablized or reduced in strength, however, in various situations. For exmaple, 
if the viscous stress is proportional only to the gas pressure, p g , i.e., r = —ap g , and a is a 
constant (see, e.g., Lightman & Eardley 1974; Coroniti 1981; Stella & Rosner 1984), then the 
disk is always stable. On the other hand, if a = ao(H/R) n (e.g., Meyer & Meyer-Hofmeister 
1983), where H is the local scale height of the disk, R is the distance from the compact object, 
and ao and n are constants, then the disk may be mildly unstable (Milsom, Chen, & Taam 
1994). The radial advection also has a stabilizing effect on the disk. This can be understood 
through the so-called S-shaped curve on the mass accretion rate and disk surface density plane, 
where the upper-branch represents the advection dominated solution. On the upper-branch, the 
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disk is thermally and viscously stable. Detailed accretion disk models with advection have been 
constructed, for example, by Abramowicz et al (1988, 1995), Kato, Honma, & Matsumoto (1988), 
Chen & Taam (1993), Narayan & Popham (1993), Narayan & Yi (1994, 1995a,b) and Chen (1995). 
The stabilizing effect of coronal dissipation on the thermal- viscous instability of accretion disks 
has been studied by Ionson & Kuperus (1984) and Svensson & Zdziarski (1994). 

In this paper, we study accretion disks with both radial advection and coronal dissipation 
included and investigate the instability behavior of such disks. In the next section we outline the 
assumptions and approximations underlying the models. The numerical results are presented in §3 
and their possible relevance to the low-frequency (~ 0.04 Hz) quasi-periodic oscillations (QPOs) 
observed in BHCs is discussed in the final section. 



2. FORMULATION 



We investigate a disk-corona system with radial advection included and focus on the structure 
and stability of the cold accretion disk instead of the hot corona above it. The system is assumed 
to be axisymmetric and non self-gravitating. The disk is optically thick and geometrically thin so 
that it can be described by the vertically integrated equations. The corona is hot and optically 
thin. Following Svensson & Zdziarski (1994), it is assumed that all the mass accretion and angular 
momentum transport take place in the cold disk. Thus in the steady state, the mass and angular 
momentum conservation equations can be written as 

M = -2irRZv r , (1) 

and 

/,£ = ^fg~\ (2) 

37T 

where M, X, v r , and v are the mass accretion rate, the surface density, the radial velocity 
and the kinematic viscosity respectively. The correction terms are / = 1 — 9ft*/(ftr 2 ) and 
g = -(2/3)(dliitt/dliiR), where r = R/R G , ft* = Sl(3R G ), and R G = 2GM/c 2 is the 
Schwarzschild radius. We assume a pseudo-Newtonian potential (Paczyhski & Wiita 1980), 



$ = -GM/(R - R G ), and the disk is rotating at the Keplerian velocity ft = y / GM/R(R - R G ) 2 . 
The viscosity is discribed with the standard a-model prescription as (Shakura & Sunyaev 1973): 

2 

v = -ac s H, (3) 

where a is a constant and c s and H are the local sound speed and the half-thickness of the disk 
respectively: 

c s = \fpfp, H = c s /ft fc . (4) 

Here p = T,/2H is the density and p is total pressure. The steady state energy conservation 
equation for the cold accretion disk is represented by the balance between the local viscous heating, 
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Q + , the local radiative cooling in the vertical direction, Q_, and the global heat transport (radial 
advection), Q a dv It is expressed as 

Q+ = Q- + Qadv (5) 

Here, irradiation heating from the corona has not been included (Svensson & Zdziarski 1994). 
This is justified by the fact that the heating of the disk by the corona only takes place in the 
surface layers down to a few Thompson optical depths (see Ross & Fabian 1993). We assume 
that a fraction tj of the total gravitational energy released is dissipated in the corona (Haardt & 
Maraschi 1991). Thus, within the disk, the viscous heating rate per unit area becomes 

Q + = (l-r,)™n 2 fg. (6) 

The local radiative cooling rate per unit area, in the optically thick radiative diffusion 
approximation, is 

4acT 4 , s 

«- = w < 7 > 

where T is the mid-plane temperature of the disk and k is the opacity assumed to be electron 
scattering, K es = 0.34. The advection cooling term is taken in a form (see, e.g. Chen & Taam 
1993): 

Y<v r p M p 

Qadv — ~p, C — Tj fyj~ S5 (8) 

R p L-kR z p 

and 

4 -3/3 din T , ^dln£ , s 

£ = -iv™j +(4 - w to (9) 

where (3 = p g /p is the ratio of the gas pressure to the total pressure, ^ = 

1 + (4 - 3/3)(7 - l)/[/3 + 12(7 - 1)(J3 - 1)], and 7 = 5/3 is the ratio of specific heats. For 

optically thick accretion disks, the equation of state is given by 

IZpT aT 4 

p= — + ~r, (10) 

fJi o 

where p is the mean molecular weight assumed to be 0.617. 



3. THERMAL EQUILIBRIA AND STABILITY 



By replacing equation (9) with a constant value of £, the surface density and mid-plane 
temperature of the disk can be solved from a set of algebraic equations for a given mass accretion 
rate and a fixed radius. This is the procedure used by Abramowicz et al (1995) to calculate the 
thermal equilibrium curves of either optically thick (see also Abramowicz, Lasota, & Xu 1986) or 
optically thin accretion disks. In the former case, the S-shaped relation between the mass accretion 
rate and the surface density is close to that calculated from the full set of slim disk equations (see 
Abramowicz et al 1988; Chen & Taam 1993). In the later case, this procedure is also confirmed by 
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Chen (1995) to be very adequate. Accordingly, we use the same approach as that of Abramowicz 
et al (1986, 1995) to examine the M(S) relation. The numerical value of £ can be estimated when 
the advection cooling is a dominant process (only then it becomes important). In that case, by 
combining equations (2)-(4), (6) and (8), one has 

£ = 44.368r- 1/2 ma- 1 #- 2 g(l - t?)" 1 ^, (11) 

where q = 1 — r _1 , m = M/M c , and M c is the critical mass accretion rate defined as 
M c = 64irGM/(cK es ). Here we have assumed 1/16 for the efficency of the conversion of rest-mass 
energy into the radiation. If we assume the disk is radiation pressure dominated then a relation 
between T and X can be derived from equations (2)-(4), and (10) as, 

T 4 = 2.464 x 10 30 r- 9 / 4 S 1 /2 rh i/2 Q ,-i/2 m -i / i/2 fl -i/2 g -3/2 j (12) 

where m is the mass of the central object expressed in unit of the solar mass. If we let the 
correction terms, /, g, and q be constants independent of radius, then X oc r -1 / 2 and T cx r -5 / 8 . 
Thus, with (3 = 0, we have, 

d In T d In X . 
dlni? dlni? v ' 

For definiteness we hereafter assume that m = 10, a = 0.1, and £ = 5.5. 

We first present the thermal equilibria of accretion disks with different tj, the fraction of 
energy dissipated in the corona (Fig. 1). For tj < 0.99, there are three branches. The lower and 
middle branches represent the local radiative cooling solutions (see Svensson & Zdziarski 1994). 
The former one is gas pressure dominated and is thermally and viscously stable. The later one 
is radiation pressure dominated and is thermally and viscously unstable. The upper-branch 
represents the advection cooling solution, it is thermally and viscously stable (see Abramowicz et al 
1988; Kato, Honma, & Matsumoto 1988; Chen & Taam 1993; Narayan & Popham 1993; Narayan 
& Yi 1994, 1995b) By increasing tj, the efficiency of the cooling (both radiative and advection) is 
effectively increased by a factor of (1 — t/) -1 . The lower and upper branches move down-right since 
Q- < Q + or Qadv < Q+ below these two lines. The middle-branch moves in the up-right direction 
since Q_ < Q + at the right side of this line. The net result is that, both the lower turning point 
and the upper turning point move to the up-right direction and \Mi — M\ \ and |Si — S2I decrease. 
Here, (Si, Mi) and (£2,-^2) denote the surface densities and the mass accretion rates at these two 
turning points respectively. Furthermore, the middle-branch becomes steeper (more perpendicular 
to the S axis) and eventually disappears for r] > 0.99. A steep middle-branch reduces the strength 
of the instability (see Milsom, Chen, & Taam 1994; Chen & Taam 1994). This result is mainly 
due to the advection, which can be illustrated through the dotted curves where no advection is 
included (Svensson & Zdziarski 1994). It is seen that the dotted lines do not turn to the right to 
form the upper-branches. The slopes of the dotted middle-branches are basically unchanged for 
various r] and they are flatter than the corresponding ones with advection included. 

The global temporal variability of accretion disks depends not only on the local instability at 
a specified radius but also on the spatial size of the disk in which it is unstable (Chen & Taam 
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1994). For the later reason, the thermal equilibria of accretion disks at different radii are shown 
in Figure 2 for two cases with tj = and tj = 0.95 respectively. For the former case, only radial 
advection but no coronal dissipation is included. It is seen that instability (the middle-branch) 
exists in wide ranges of both disk radii and mass accretion rates. This results in a global, burst-like 
time-evolution (see Taam & Lin 1984; Lasota & Pelat 1991). The inclusion of coronal dissipation 
apparantly weakens the instability of the disk in three aspects (in comparision to tj = 0). First, 
the spatial size where the disk is locally unstable is reduced to r < 100 from r > 500. Second, the 
range of mass accretion rates at which the disk is unstable is greatly reduced to a neighborhood of 
~ 0.1 — 1M C . Third, the strength of the local instability is also reduced (steeper slope). Thus the 
global temporal variability of the disk can be very mild or totally disappears for tj > 0.95. 

4. DISCUSSION 

We have demonstrated that accretion disks with coronae may have only a very mild 
thermal- viscous instability in comparision to models without a corona. This is due to the 
stabilization effects of coronal dissipation and radial heat transport within the cold disk. We 
suggest that this result is applicable to the low-frequency QPOs (~ 0.04 Hz) observed in the low 
states of black hole candidates Cyg X-l (Angelini k White 1992; Vikhlinin et al 1992a, 1994; 
Kouveliotou et al 1992a) and GRO J0422+32 (Vikhlinin et al 1992b, Kouveliotou et al 1992b, 
Pietsch et al 1993). 

It is well-known that accretion disks which possess an S- Shaped thermal equilibrium curve 
in the mass accretion rate and surface density plane (S,M) may experience a limit-cycle time 
evolution. This type of thermal- viscous instabilities has been widely studied by using nonlinear 
time-dependent calculations for accretion disks surrounding either white dwarfs or neutron stars 
and black holes. Even though the physics of the origin of these instabilities is different, their 
evolution behavior is quite similar (see Cannizzo 1993 for a review for dwarf novae and Taam & 
Lin 1994 and Lasota & Pelat 1991 for neutron stars). For example, for burst-like evolutions, it 
is understood that the recurrence time of bursts corresponds to the viscous diffusion time of the 
regions where the disk is unstable. On the other hand, the rising time of the burst corresponds 
to the thermal time scale and which is shorter than the viscous time scale. Thus, by reducing 
the radial width of the unstable disk the evolution may become oscillation-like once the thermal 
and viscous time scales become comparable. This is exactly the case as was shown by Chen 
& Taam (1994) based on a modified viscosity prescription. It was also shown there that the 
amplitude of the oscillation becomes smaller when the strength of the instability is weaker, which 
is represented by a steeper slope of the middle-branch of the thermal equilibrium curve. A type 
of small amplitude, periodic (~ 0.04 Hz), oscillations has been obtained for both neutron star 
and black hole accretion disks and which was suggested to be relevant to the low-frequency QPOs 
observed (Chen & Taam 1994). For the present model, similar global time-dependent calculations 
are needed to constrain the disk-corona parameters, such as a, tj, m, and M, quantitatively. 
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Nevertheless, the nature of the instability implies a time scale of R 2 jv ~ 10 5 a 1 (R/ H) 2 r 3 / 2 s, 
which is in the order of tens of seconds for r ~ 100 and H/R ~ 0.01 — 0.1. 

We have adopted the concept of coronal dissipation popularly applied in the recent two-phase 
disk-corona models (Haardt & Maraschi 1991). Ionson & Kuperus (1984) have used a slightly 
different approach. In their assumption, the disk consists of an internal and an external viscosities; 
only the internal one generates heating in the disk, while both of them transport the angular 
momentum. Since the physics of tj is not well understood, we leave it as a free parameter. We 
confirm the result of Ionson & Kuperus (1984) and Svensson & Zdziarski (1994) that coronal 
dissipation has a stabilizing effect on the thermal- viscous instability of accretion disks by raising 
the mass accretion rate threshold above which the disk is unstable. 

We have not included the irradiation heating from the corona on the disk. The irradiation 
heating can dominate the viscous heating in the disk if it is deposited throughout the cold disk 
(see, e.g., Nakamura & Osaki 1993; Kusunose & Mineshige 1994). However, since the deposition 
and heating of the radiation from the corona only takes place in the surface layers down to a few 
Thompson depths (see, e.g., Ross & Fabian 1993), irradiation heating will probably not affect 
the cold optically thick disk very much (Svensson & Zdziarski 1994). Furthermore, it has been 
also known that the effect of irradiation is unimportant in the deep interior of the disk so long as 
the irradiated flux is much less than the product of the viscously generated flux and the optical 
depth of the disk (Lyutyi & Sunyaev 1976; Taam & Meszaros 1987). In the disk-corona system, 
the irradiation (from corona) and viscously generated (from disk) fluxes scale to tj/2 and 1 — tj 
respectively. Thus the criterion for unimportant irradiation reads kT, — rf). This is easily 

satisfied even for tj = 0.99 due to a large surface density, X ~ 10 5 — 10 6 . 

The structure and stability of the corona is beyond the scope of the present paper. In fact, 
the corona should not be examined separately from the disk-corona system. In a self-consistent 
disk-corona model, one may need to consider the disk-corona system together and should probably 
address problems such as the radiation transfer and its two-dimension (radial- vertical) nature. 
Futhermore, the non-Keplerian effect and the transonic nature of accretion flows near the vicinity 
of the black hole will also play important roles in determining the global structure of the disk. 
These will be the subjects of our future investigations. 

The author would like to thank Profs. Marek Abramowicz and Ronald Taam for discusions and 
also an anonymous referee for useful suggestions. 
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log (S) 

Fig. 1. — The S-shaped curves (solid lines) of M(S) relation for various r], the fraction of energy 
dissipated in the corona. Here r = 10, a = 0.1, and m = 10. Note that for tj > 0.99, the entire 
disk is locally thermally and viscously stable. The dotted lines represent the solutions without 
advection. 
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Fig. 2. — The S-shaped curves of M(S) relation at different radii. Note that for r] = 0.95, the 
unstable behavior of the disk is restricted to a relatively narrow spatial region and is present in 
a small range of mass accretion rates. The dotted lines in b) represent the solutions without 
advection. 



